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ABSTRACT 

A 30-day contiguous photometric run with the MOST satelhte on the WN5- 
6b star WRllO (HD 165688) reveals a fundamental periodicity of P = 4.08 ± 
0.55 days along with a number of harmonics at periods P/n, with n 2,3,4,5 
and 6, and a few other possible stray periodicities and/or stochastic variability 
on timescales longer than about a day. Spectroscopic RV studies fail to reveal 
any plausible companion with a period in this range. Therefore, we conjecture 
that the observed light-curve cusps of amplitude ~ 0.01 mag that recur at a 
4.08 day timescale may arise in the inner parts, or at the base of, a corotating 
interaction region (CIR) seen in emission as it rotates around with the star at 
constant angular velocity. The hard X-ray component seen in WRllO could 
then be a result of a high velocity component of the CIR shock interacting with 
the ambient wind at several stellar radii. Given that most hot, luminous stars 
showing CIRs have two CIR arms, it is possible that either the fundamental 
period is 8. 2d or, more likely in the case of WRllO, there is indeed a second 
weaker CIR arm for P = 4.08d, that occurs ~two thirds of a rotation period 
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after the main CIR. If this interpretation is correct, WRllO therefore joins the 
ranks with three other single WR stars, all WN, with confirmed CIR rotation 
periods (WR 1, WR6, and WR 134), albeit with WR 110 having by far the lowest 
amplitude photometric modulation. This illustrates the power of being able to 
secure intense, continuous high-precision photometry from space-based platforms 
such as MOST. It also opens the door to revealing low-amplitude photometric 
variations in other WN stars, where previous attempts have failed. If all WN 
stars have CIRs at some level, this could be important for revealing sources of 
magnetism or pulsation in addition to rotation periods. 

Subject headings: stars: individual (WRllO) — stars: winds, outflows — stars: 
rotation — stars: Wolf-Rayet 



Introduction 



While most adequatel y observed O stars exhibit pairs of opposing Corotating Inter- 
action Regions (CIRs; see iMuUan Ill984l ). manifested by propag ating Discrete Abso rption 



Components in their resonant P Cygni absorption components (iKaper et al. Ill999l ). little 



such evidence is found for CIRs in 0-star descendants, the Wolf-Rayet (WR) stars. If found 
in WR stars, CIRs would give precious information on their rotation, which would be highly 
relevant in the context of rapi dly rotating WR stars as prime candidates for slow Gamma- Ray 



Bursts fISt-Louis et al. Il2009[ ). 



Recently ISt-Louis et al. I ( l2009l ) (see also IChene fc St- Louis I (120101 )) have examined a 
northern sample of 25 WR stars for the presence of CIRs via large-scale variations exhibited 
on broad, strong emission lines. Some 20% of these stars show clear signatures of such effects, 
with WR 110 situated at the limit between those stars that show CIRs and those that do not. 
In addition, some (and possibly all, given the limited number of spectra) stars in the sample 
unsurprisingly show small-scale wind-clump variations. Such variations manifest themselves 
as stochastic pro pagating small-scale sp ectral variations on time series of virtually all WR 
line spectra (e.g. iLepine &: Moffat 1119991 ). 



WRllO is a moderately bright {y = 10.30,6 — v = 0.75), reddened {A^ = 3.83 mag) 
southern [RA(2000) = 18:07:56.96, DEC(2000) = -19:23:56.8] Galactic (/ = 10.80°, b = 



^Based on data from the MOST satellite, a Canadian Space Agency mission, jointly operated by Dynacon 
Inc., the University of Toronto Institute for Aerospace Studies and the University of British Columbia, with 
the assistance of the University of Vienna. 
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+0.3 9°) Wolf-Rayet (WR) star in the mid-nitrogen sequence with type WN5-6b (jSmith et ah 



19961 ). The "b" suffix indicates broad hues for its subtype and thus low H-a bundance by de- 



fault . Located at a distance of 1.28 kpc, it is 9 pc above the Galactic plane (Ivan der Hucht 



200ll ). in the normal range for most extreme population I stars. Although to date it lacks 



any obvious binary signature, a comprehensive radial velocity (RV) search has not yet been 
carried out. However, WRllO does have relatively high X-ray flux, with a significant hard 
component above 3 keV compa red to its lower-energy 0.5 keV emission, as seen in most 
WN stars (jSkinner et al. II2002I ). The lower-energy component likely arises due to small- 
scale turbulent shocks in the wind (assuming thermalization of matter colliding at velocity 
v: kT l/2mpv'^), with veloci ty dispersion ~ 10^ km s~^, as seen i n virtually all WR-star 



winds ([Lepine &: Moffat I Il999l ) and other luminous hot stars (e.g. lEversberg et al. I Il998 



Lepine &: Moffat I l2008l ). The former hard component implies velocities of ~10^ km s 



which are more difficult to explain, without invoking an additional mechanism to provide 
such high speeds (e.g. accretion onto a low-mass companion or colliding winds with a mas- 
sive companion]. H owever, WRllO emits normal thermal radio emission from its wind 
(jSkinner et al. 1120021 ). which lends little support for its binary nature. As we show in this 
investigation, it is possible that CIRs are present in the wind of WRllO, and that WRllO's 
hard X-rays may arise in the shocks produced by its CIRs at velocities of ~10^ km s^^. 



2. Observations 

The optical observations of this work consist of both high-precision, nearly continuous 
photometry over a month {Tph = 30 days) from space using the Microvariability and Oscil- 
lations of STars (MOST) microsatellite, and ground-based optical spectroscopy from various 
observatories over 8 years. 



2.1. Photometry 



MOST (IMatthews et al. II 19991 : 1 Walker et al. 1120031 1 contains a 15-cm Rumak-Maksutov 
telescope imaging onto a CCD detector via a custom optical broadband (350-750 nm) filter. 
From its polar Sun-synchronous orbit of altitude 820 km and period 101 min, MOST has 
a continuous viewing zone (CVZ) about 54° wide within which it can monitor target fields 
for up to two months without interruption. Targets brighter than V ~ 6 mag are observed 
in Fabry imaging mode, while fainter targets (like WRllO) are observed in direct imaging 
mode (3" x 3" pixels in a ~ 1° field), similar to standard CCD photometry with a ground- 
based instrument. The photometry is non-differential, but given the orbit, thermal and 
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design characteristics of MOST, experience has shown that it is a very photometrically stable 
platform even over long timescales (with repeatability of the mean instrumental flux from a 
no n- variable target of the brightness of WRllO to within about 1 mmag over a month). 

Due to the failure of the tracking CCD, the science exposures on the MOST science CCD 
now take place at the same time as the guide-star exposures for satellite pointing. In the 
case of the WRllO observations, the guide-star exposure time (and hence also the science 
target exposure time) was 0.3 sec. To build up sufficient signal-to-noise ratio (S/N), the 
science exposures were co-added on board the satellite in "stacks" of 100 exposures. Stacks 
(each 30 sec of total integration) were downloaded from the satellite consecutively (with no 
dead time between stacks), giving a sampling rate of about twice per minute. 

The WR 110 field lies outside the MOST science CCD CVZ, so it could not be monitored 
(strictly speaking) continuously. This field alternated with another MOST primary science 
target (Altair) during each MOST satellite orbit. Typically, we monitored WRllO for about 
40-55 % of each MOST orbit (i.e. about 40-55 min of each 101 min) during the low-straylight 
portion of the orbit, with ~ 35 000 individual 30-second data-points between 2007 June 17 
and July 16 (HJD 2 454 267 - 297). 

Figure [1] shows the entire MOST light curve, binned for each MOST orbit. Although 
significant variations are seen on timescales longer than ~a day, few variations are seen on 
shorter timescales and binning allows one to increase the photometric point-to-point rms 
precision from ~3 mmag before binning to ~0.5 mmag after binning. Interpretation of the 
observed variations is given below in the next section. 

2.2. Spectroscopy 

Nine CCD-spectroscopic runs were carried out at five telescopes over eight years, in- 
cluding one parallel run with the MOST observations. However, this run was plagued by 
bad weather and instrumental problems, resulting in limited time coverage. A journal of 
observations is given in Table [TJ Generally, the spectra cover timescales ranging from ~10 
minutes to several years and range from near the Balmer limit to HelA5876, although not 
the same for each telescope. A SNR~100 was obtained per extracted pixel, which ranged 
from 0.5 to 1.6 A, with resolution of 2-3 pixels. The bias subtraction, fiat-fielding, spectrum 
extraction, sky subtraction and wavelength calibration of all spectra were executed in the 
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usual way using IRAFg software. 

Special care was taken for the normalization of the spectra. First, a mean was made for 
each run. Then each spectrum of a run was divided by the run mean and the ratio fitted 
with a low order Legendre polynomial (between 4*^^ and 8^^ order). The original individual 
spectrum was divided by this fit, and was therefore normalized to the run mean. When 
this procedure was done for each run, the run means were then put at the same level by 
using the same procedure as described above. This allowed us to put all individual spectra 
at the same level. Then, we combined all run means into a global, high quality mean, 
which was then fitted in selected pseudo-continuum regions, i.e. wavelength regions where 
large emission-lines do not dominate. These regions are shown in Figure |2j Finally, the 
fitted continuum function is applied to rectify each individual spectrum. The error on the 
continuum normalization - measured as the standard deviation of individual spectra around 
the continuum function - is typically of 0.5%. Mean spectra from each telescope are shown 
in Figure [2l 



3. Analysis 
3.1. Light Curve 

We start by examining the Fourier spectrum for the MOST photometry (Figure |3]). 
Here the largest amplitude occurs at a frequency of z/q = 0.245 ± 0.033 cycles day~^ - with 
error equal to l/Tph, where Tph is the total duration of the light curve - (P = 4.08 ± 0.55 
days), with harmonics of this frequency at, or close to, n x uq, with n ~ 2, 3, 4, 5 and 6. A 
few other peaks do occur at low (<20% of the dominant peak) relative amphtudes, notably 
at periods longer than ~a day, but with no obvious order or interrelationship. None are 
due to sidelobes of the primary maximum, due to the uniform, dense data sampling. No 
significant peaks beyond the above harmonics were detected with periods shorter than ~a 
day. We assume that this implies a single dominating (non-sinusoidal) periodicity with P 
= 4.08 days, plus mainly stochastic variability on timescales longer than ~a day (possibly 
related to wind clump activity, which survives for at least ~8 hours in typical WR winds: 



Lepine fc Moffat Ill999l ). Examination of the light curve in Figure [T] shows relatively sharp 
cusps occurring on this time scale (i.e. P = 4.08 days, starting at to = HJD 2 454 270.27) 
in all seven cases but one (the third case), although even then there is a rise in brightness, 



^IRAF is distributed by the National Optical Astronomy Observatories (NOAO), which is operated by the 
Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the 
National Science Foundation (NSF). 
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though not cusp-hke. The cusps have typical widths at the base of ~a day. Such cusps are 
clearly non sinusoidal, requiring a series of successive harmonics to reproduce. Figure [3] also 
shows a time-frequency plot with straddled 8-day (double the best period) sampling. This 
plot reveals that the 4.08-day period stands out at all times except near the middle and at 
both ends (where edge effects come into play) of the photometric time series. The former is 
compatible with the apparently suppressed third cusp, as noted above. 

Figure m shows the phase-binned light curve of the orbital-mean MOST data from Fig- 
ure dl based on the period of 4.08 days and zero point to (nearly coinciding with the first 
observed cusp in Figure [1]). This light curve is quite noisy compared to the appearance 
of the non phase-binned light curve in Figure [TJ Further binning to 0.05 in phase reveals 
a dominating peak at phase zero, with a second lower peak at phase ~0.65. Both peaks 
can be traced back in Figure [T] often to sharper, clearer cusps at similar repeated intervals. 
These cusps become somewhat smeared, although remain still visible in Figure HI probably 
due to dispersion in shape and timing around some average, combined with stochastic wind 
variations. 



3.2. Spectral Line Variations 

3.2.1. Short-term 

Only one spectroscopic run (CTIO 2004) is of sufficient quality and length to allow a 
search for variability on timescales less than a day. Figure O shows greyscale time-plots of the 
four successive CTIO nights for differences from the mean of each spectrum, for the principal 
WR spectral emission lines of HellA5411 and ClvAA5802/12 + HelA5876. From these plots, 
one sees evidence for the usual stochastic, emission subpeak s, propagating outward s from 



line center, as seen in all previously well-studied WR spectra (ILepine fc Moffat Ill999l ). Such 
perturbations are believed to arise in turbulent wind clumps, with no apparent periodicity 
on timescales of a night (~0.3 day) or shorter. 

The simplest, most effective way to examine spectral line variations is via time depen- 
dence of line moments from order zero through four. The zeroth moment gives the line 
strength or equivalent width (EW) when normalized to the continuum. The first moment 
yields the global line position or radial velocity (RV). The second moment can be used to cal- 
culate the line width a la gaussian. Line asymmetries can best be estimated by (normalized) 
third moments, known as skewness. Finally, normalized fourth moments yield information 
on how the profile differs in global shape from a gaussian. 



Figure [6] shows a time plot of normalized EW, RV, skewness and kurtosis of the CTIO 
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observations for HellA5411. The EWs were calculated by integrating the function (1 — 
Fx), where Fx is the rectified line flux between AA=5348-5460 A. In order to highlight 
the variability, we prefer to divide the EWs by the mean value of all EWs. The radial 
velocities (RVs) were obtained by cross correlating individual spectra in the wavelength 
interval AA=5140-5980 A with the mean spectra. (Line bisectors were also explored to obtain 
RVs, but they revealed larger scatter yet similar trends as the cross-correlation technique.) 
The other moments of the HellA5411 line are calculated in the wavelength interval AA=5348- 
5460 A. The n^^ central moment is defined as follows: 

/i„ = S,(A,-A)"/,/S,/,- (1) 

where 

A = S,A,.J^./S,/„ (2) 

with Ij being the intensity of the line and Xj the wavelength. The skewness is jJ.'i/ ^^2'^ and the 
kurtosis is fJ^^/fJ^"^- Variability is seen in all of these quantities, although with no significant 
periodicities. Any longer-term periodic variations of these data will be examined together 
with the remaining spectra in the next subsection. 



3.2.2. Long-term 



Figure [7] shows time plots of EW, RV, skewness and kurtosis (measured as above) for the 
strongest, least blended lines: (a) over the full 6 years for HellA4686 and (b) T^p = 8 years for 
HellA5411. In the case of the latter from CTIO, we have taken nighly averages, so that the 
numerous original CTIO data will not dominate the period analysis. The moderately strong 
lines of ClvA5802/12 and HelA5876 are blended and seriously affected by telluric lines, so 
they were not analyzed. 

While these four parameters all show variations for both lines, those of largest amplitude 
are in RV and skewness for HellA4686. These two quantities appear to vary in an anti- 
correlated way (correlation amplitude, r = -0.89), as expected if the line varies in shape; e.g. 
an excess on the red side of the line will cause a small positive RV shift and a negative change 
in skewness. A similar trend (r = -0.86), is also measured for these two parameters in the 
HellA5411 line, despite the apparently lower amplitude of the variations. Presumably, these 
are variations in profile shape and do not reflect tr ue RV motions of t he star. Furthermore, 
a search using phase-d ispersion modulation (PPM. IStellingwerf Ill978l ) and Fourier methods 
( iLenz fc Breger 1 12005[ ) was carried out for periodicities in the RVs in the range expected 
from the MOST variations This was done in the frequency range u = 0.2037 - 0.2885 cycles 
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day~^, covering the whole photometric Fourier peak centered at P = 4.08 days, and with 
frequency step size 2 x 10^^ cycles day~^, i.e. ~ (10 x 2Tsp)~^. The search failed to reveal 
any significant periodicity. Alternatively, a phase RV plot of each line phased with the best 
period in this interval produces a scatter diagram with force-fitted sine-wave RV amplitude 
of K ~14, 10 and 28 km s~^ for all RV data (nightly means for CTIO, with typical errors 
of ~10 km s-i) for HellA5411, the 335 CTIO data alone for HellA5411, and all RV data 
for HellA4686, respectively. The lower the scatter, the lower the K value. Combined with 
the above information on anti-correlation of RV and skewness, we deduce that if WRllO 
is a binary system, it is not evident from the data presented here. On the other hand, the 
lack of a clear indication of a 4.08 day period in any of the spectroscopic data does not 
necessarily mean that such a period is not present. We suspect that it will take a much 
larger, higher- quality spectroscopic sample to reach any firm c onclusions on any p eriodicity 
in spectral behavior of WRllO, as was the case with WR134 (IMorel et al. Ill999l ). 



4. CIR model 

Plausible scenarios to account for the 4.08 day periodicity include binarity, pulsations 
and CIRs. Binarity is unlikely, due mainly to a lack of coherent RV variations that phase with 
the primary photometric period. Pulsations are unlikely due to the highly non-sinusoidal 
form of the variations at P = 4.08 days. Hence, the most likely scenario to explain the 
low-amplitude (~1 %) 4.08 day periodicity and cusp-like behavior in the MOST light curve, 
along with mostly unrelated spectral variations, appears to be large-scale, emitting over- 
density structures rotating with the wind. (One also has the impression of cusp-like minima 
in the light curve; however, these could be a result of the small space between two emitting 
cusps occurring per rotation - see below.) This would be similar to what is seen in a few 
other WR stars (WR6, WR134, WRl - see below). Corotating interact ion regions (CIRs; or 



their s ource at the base of the wind) as proposed for O star winds by ICranmer &: Owocki 



( 119961 ) provide an attractive scenario to explain our obser yations. The Cranra er & Owocki 



model is designed for O stars, explicitly assuming CAK ( ICastor et al. Ill975l ) parameters, 
which may not apply directly to H-free WR stars. Nevertheless, the idea of a CIR driver at 
the base of the wind in WR stars is a qualitatively plausible scenario, analogous to 0-star 
winds. We now apply this idea to WRllO in the following way. 

First, in the context of the Cranmer & Owocki model, we assume that the CIRs are 
continuum-emitting thermal sources, heated by the associated shock action created by a 
hot spot on the underlying rotating star interacting with the ambient wind. Most of the 
CIR emission must arise close to the star, where the ambient wind density is highest, as 
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inspired by the simulations of ICranmer &: Owocki I (119961 ) - see a simulation in their Fig. 5. 
Alternatively, the emission could arise in a hot spot at or close to the stellar surface, that 
produces the CIR. To simplify our calculations, we assume a point source located at some 
radius Rg = 7-R* from the cen ter of the star, with 7 > 1 a nd R^, being the stellar core radius. 
For simplicity (and following ICranmer fc Owocki 1 119961 ) we assume the CIR to initiate at 
the stellar rotation equator. 

In this scenario, the changes in the light-curve are caused by an assumed CIR-associated 
point source attenuated by the wind and seen at different angles with the line of sight as 
it rotates on the near side of t h e WR star. In this study, we use a simple WR wind model 
derived by iLamontagne et al. I ( 119961 ) in the context of wind eclipses for WR + O systems. 
Hence, the change in magnitudes is defined as : 



Am = constant — 2.5 log^^Q (/pv_R + /^e 



(3) 



with arbitrary constant and where Iwr is the intensity of the WR star and Ig is the intensity 
of the hot spot. Here the total opacity r between the source and the observer, passing 
through the WR wind is : 



k / diz/Rg) ir/Rgf{l-RjrY 



-1 



(4) 



where [r/Rg 



(cosi cos 27r0)^ + sin^27r0 + {z/Rg 



— sinz) cos27r0, i is the 



inclination of the rotation axis (assumed to coincide with the CIR rotation axis) relative 
to t he observer, and (/) is the rotation phase. We allow for a WR wind with a /3 law (like 



m iLamontagne et al. Ill996l ). in which the actual /3 value will be fitted. Since the CIR 
source is much closer to the inner WR wind than the O star is in the (non photospheric 
eclipsing) binary cas e, we expect that = w ill not be a good choice, as it wa s in the case 
of WR + O binaries ( ILamontagne et al. Ill996[ ). Following ISkinner et al. I (120021 ) . we adopt a 
radio-based thermal mass-loss rate for the WR star of M = 1.6 10~^ Mq yr~^, after dividing 
by a factor three to allow for wind clumping, a wind terminal speed o f f no =2100 km s~^, 
and a WR core radius i?* = 4 Rq. We find an opacity constant as in ILamontagne et al. 
(119961 ) k = aMae/[4:nmpVooRs] = O.28/7, where a = 0.5 is the number of free electrons per 



He nucleus in a highly ionized wind. In our case, we are primarily interested in the phases 
centered at phase 0.50 (i.e. when the source is between the star and the observer) since the 
CIR point source will not be normally visible as it rotates on the back side of the WR star. 

We calculated Equation [3] numerically, centered at phase 0.5 for half a rotation period, 
as a function of the parameters (3, Ig/IwR, h and 7. We consider two cases: (a) an optically 
thin point source, whose emission is isotropic in direction, and (b) an optically thick source 
whose emission is maximum when seen perpendicularly to the stellar surface, dropping to 
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zero when seen parallel to the surface (i.e. like pancakes on or near the stellar surface). In 
case (a) light-curve excesses originate from the hot spot and in case (b), we have to multiply 
the net spot intensity (J^) from Equation [3] by the projection factor (sini) cos[27r(0 — 0.5)]. 

Taking the second, the fourth and the fifth cusps in Figure [1] as the cleanest, most 
representative for m of the hotspot light curve, we fit the above parameters using the IDL 
routine mpfit.pro (IMarkwardt II2009I ). The fit does not converge when the inclination i is set 
as a free parameter. Therefore, we fix it to different angles in the range i = 90° — 40° with 
increments of 10° before fitting all the other parameters. A summary of fitted parameters 
using different inclination angles is given in Table [2j Note that we do not include our results 
for an inclination i lower than 40°, since the fit does not converge when using such small 
angles. The values of Ig/IwR tend to be slightly smaller in case (b) than in case (a), due to 
the projection factor in case (b) that requires less extended wind opacity to bring about the 
observed cuspy shape. Indeed, it is the varying opacity through the WR wind as the WR 
star rotates, which leads to the cuspy shape; the closer the source is to the stellar surface (i.e. 
7 closer to unity), the smaller the (3 value (and less extended the wind) necessary to give 
the cuspy shape. The best solutions for the thin and the thick cases at different inclination 
angles are plotted in the left panel of Figure [3 We also repeated the fits for the thin and 
thick cases at the different inclination angles using only the second cusp, in order to verify 
how our results would vary depending on this choice. The best solutions are plotted in the 
right panel of Figure [8] and also summarized in Table [2l The difference between the results 
using one or three cusps may come from an epoch-dependency of the variations, smaller cusps 
present in only some of the considered cusps and/or a change in the shape and amplitudes 
between the different cusps due to superposed stochastic variability. 

All the fits give a Xred close to unity. It is hence, at first, not possible to exclude any 
result only based on the quality of the fit. As seen in Figure [H], changing the inclination does 
not change the shape of the curve as much as the amplitude and requires higher values of /3, 
7 and Is/Iwr (this may explain why the fit does not converge when the inclination i is set 
as a free parameter). However, in many cases Is/Iwr > 1? which means that the intensity 
of the spot is greater than the intensity of the whole star. Since such a solution is not likely, 
we may exclude the optically thin case, and set i > 70°. The remaining solutions imply 
(3 ~ 11.5 — 12.5. These value would give (3R^ ~ 5Oi?0, which is comparable to the value 
obtained by iLepine fc Moffat I (119991 ) for stars of similar spectral type. 

Although the overall fit is quite good, the model does not reproduce well the peaked 
central cusps. Tweaking the model to do this may require assumptions that go beyond 
the scope of this study. For example, it might b e possi ble that the expected density- deficit 
channel parallel to the CIR ICranmer &: Owocki I (119961 ) could affect the light-curve peaks. 
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Discussion and Conclusions 



Assuming that CIRs do in fact account for the variabihty observed in the MOST hght 
curve of WRllO, we can now provide a plausi ble explanation for the har d X-rays coming 
from WRllO. According to the models of CIR (jCranmer fc Owocki Ill996l ). shock velocities 



in the radial direction can reach over 10^ km s~^ at several stellar radii (beyond where 
the optical hot spot is seen), although highly variable from one m odel to another. This is 
several times larger than the velocity dispersion (~ 10^ km s~^: iLepine fc Moffat 1 119991 ) 
in clumps that is likely associated with the low-energy X-rays (kT 0.5 keV). Scaling up 
proportional to the square of the velocity, one can easily account for the observed hard 
X-rays using CIRs without accretion onto a compact or low-mass binary co mpanion. Note 
that the wind terminal velo city of WR 110, altho ugh high (2100 km s~^ from lSkinner et al 



2OO2I . or 2300 km s~^ from lHamann et al. II2006I ). is not the source per se of the X-rays; it 
is the velocity dispersion or the radial CIR velocity component which counts for producing 
the low- versus high-energy X-ray shocks, respectively. 

Our modeling of the photometric light curve leads us to conclude that WRllO joins 
the ranks of three other WR stars, all WN, that exhibit CIRs with known rotation periods 
of the underlying sta r: WR6, WN4, P = 3.76 days; WR134, WN6, P = 2.3 days; WRl, 
WN4, P = 16.9 days jMorel et al. Ill997[ll999l : lchene fc St-Louis Iboiol . respectively). While 
the periodicity in WR6 was easy to find because of the high photometric amplitude (~ 0.1 
mag), WR134 was more difficult with ~ 0.03 mag plus significant noise from other sources, 
and despite its large amplitude ~ 0.1 mag, WRl has a long period and complex light curve. 
WR 110 has an even lower amplitude, with significant noise, making the availability of MOST 
for intense, long-term monitoring indispensable for the discovery of its CIRs. In all 4 cases, 
the cuspy nature of the often multiple light-curve bumps, is evident and may thus be due 
to the same explanation: continuum-emitting CIRs that emit mainly just outside the inner 
opaque WR wind, where the intensity falls off rapidly at larger radii. In mos t cases, there 



are 2 cusps, implying 2 CIRs. This is the usual number found in 0-star winds (IKaper et al. 
1999h . 



None of these rotation rates appear to be rapid enough to be directly associate d with 



Gamma-ray Bursts, assuming rotational coupling between the surface and the interior (ISt-Louis et al. 
2009h . 



Despite its very small amplitude and complex pattern, WRllO has possibly revealed its 
CIR rotation properties due to the power of the MOST satellite, i.e. its ability to observe 
continuously for a month at high precision. This may open up the door to revealing other 
subtle cases, giving us a better idea of the properties of WR stars and their strong winds. 
It may even have given us a lead to probing the mysterious origin of the CIR phenomenon 
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rooted at the stellar surface, e.g. as a phenomenon related to magnetism or pulsations. 
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Fig. 1.— Upper. MO^T light curve of WRllO in 101-minute (MOST orbit) bins. Tlie 
vertical lines indicate intervals of period 4.08 d starting at to — HJD 2454270.27. Lower. 
Simultaneous MOST light curve of HD 312657, V = 10.8, Sp = A2, reduced in the same 
way as WRllO, with arbitrarily displaced magnitude zero-point. This star appears to be 
stochastically variable at the 0.005 mag peak-to- valley level and slightly larger point-to-point 
rms scatter of ~1 mmag, compared to ~0.5 mmag for WRllO. 
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Fig. 2. — Mean spectra from each run as in Table [H with some key hnes identified. The 
vertical scale refers to the lowest spectrum, with the remaining spectra vertically displaced 
for clarity. The insert shows a vertical compression of the strong HellA4686-line region. The 
horizontal bars show positions used to calculate the continuum. 
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Fig. 3. — Upper panel: Fourier amplitude spectrum of the binned MOST light curve in 
Figure [H with insert showing expanded scales. The highest peak refers to the adopted 
fundamental frequency corresponding to a period of P = 4.08 days. Harmonics are indicated 
at frequencies at multiples of 2 through 6 times fundamental (marked with downward arrows). 
Lower panel: Time-frequency Fourier plot with 8-day running windows in time. 
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Fig. 4. — Phased light curve from Figured! using P = 4.08 days and to = HJD 2454270.27. 
The large points show means in phase bins of 0.05, with 2-a error bars based on the standard 
error of the mean in each bin. 
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Fig. 5. — Upper 4 panels: Grey-scale plots of the difference from the mean of the CTIO 
spectra over four successive nights (shown in the bottom panel), (a) HellA5411, (b) 
CIVIVA5802/12 + HeiA5876. 
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Fig. 6. — Normalized equivalent width (EW), line radial velocity (RV), skewness and kurtosis 
versus time for the CTIO data of HellA5411. 
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Fig. 7. — As in Figure [6]for all the spectroscopic runs for (a) HellA4686 and (b) HellA5411, 
with EW relative to the global mean (excluding the OMM data, due to detector zero-point 
intensity problems). The CTIO data for HellA5411 are shown in nightly means. 
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Fig. 8. — LeftThe data points refer to the second (black diamonds), the forth (blue triangles) 
and the fifth (red squares) cusps indicated by vertical lines in Figure [H We assume that 
these cusps begin and end when the light curve reaches a minimum, i.e. when there is 
no enhancement due to the passage of a CIR. To avoid contamination from previous and 
subsequent features, three cusps are selected in the time intervals HJD 2 454 273.725 - 
275.000, HJD 2 454 282.240 - 283.450 and HJD 2 454 285.770 - 286.720, respectively. To 
facilitate the fit, the data points are shifted vertically by an arbitrary constant, so the 
minimum is fixed at Am=0. The curves indicate the best fits at different inclination angles 
(see Table [2]) for an optically thin point source, using black dashed lines, and an optically 
thick source, using blue dash-dotted lines (see text). Note that these curves are superposed 
and hard to distinguish Right: Same as the left panel, but only using the second cusp. 
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Table 1: Journal of spectrosocpic observations. 



Telescope 


Dates 




No. spectra 


Range (A) 


A/pi 


Leoncito 2.2m 


2001 Apr 8 - 


10 


3 


3900-5500 


1.63 




2002 Apr 20 


- 23 


4 


?) 


)) 




2003 Oct 5 




1 


?) 


5) 




2004 May 20 


-31 


8 


?) 


5) 




2004 Jun 13 - 


- 16 


4 


5) 


5) 


CFHT 3.6m 


2003 Aug 9 - 


13 


17 


4300-6000 


0.51 


CTIO 1.5m 


2004 Jun 29 - 


- Jul 2 


335 


5300-6200 


0.77 


OMM 1.6m 


2007 Jul 13 - 


17 


5 


4400-6100 


0.66 


Leoncito 0.6m 


2009 Mar 12 


- Apr 16 


15 


3700-5300 


0.97 
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Table 2: CIR model fit parameters. 
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9 fiSI 


n 47fi 


1 4'^n 






fin 

OU 


1 1 7fi 


9 71 fi 


n ^QA 


1419 






50 


12.08 


2.778 


0.847 


1.386 






40 


12.56 


2.863 


1.488 


1.350 


Thin 


1 


90 


14.19 


2.895 


2.595 


1.198 






80 


14.20 


2.899 


2.685 


1.194 






70 


14.30 


2.920 


2.993 


1.180 






60 


14.51 


2.960 


3.666 


1.155 






50 


14.81 


3.017 


5.151 


1.119 






40 


15.32 


3.102 


9.095 


1.068 


Thick 


1 


90 


12.48 


2.744 


0.779 
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80 


12.51 
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0.813 
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12.62 


2.773 


0.929 
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2.808 


1.187 


1.121 






50 


13.12 


2.870 


1.758 


1.094 






40 
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2.956 


3.294 


1.055 



